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Very recently the BaBar collaboration has put forward a claim that the X(3872) is not a 1 ++ 
resonance, as most of the phenomenological work on the subject was relying on, but rather a 
2 ^ one. We examine the consequences of this quantum number assignment for the solution 
of the X(3872) puzzle. The molecular interpretation appears less likely, and the conventional 
charmonium interpretation should be reconsidered. There are several well-known difficulties with 
this interpretation, to which we add a new one: the production cross section at CDF is predicted 
to be much smaller than that observed. We also confirm, using a relativistic string model, the 
conclusion from potential models that the mass of the state is not consistent with expectations. 
In the tetraquark interpretation the 2 h assignment implies a rich spectrum of partner states, 
although the V(3872) may be among the few which are narrow enough to be observable. 
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I. INTRODUCTION 



In a very recent paper the BaBar collaboration claims [1] that the quantum numbers of the A(3872) are not 
1 ++ , as had generally been accepted, but 2 '". Whereas an early conference paper by Belle [2] favored 1 ++ 
quantum numbers, a later analysis by CDF with higher statistics [3] allowed both 1 ++ and 2 h . The new 
conclusion from BaBar is based on the distribution of the 37r mass in the J/ipto decay near threshold, which 
strongly favors a P-wave rather than S-wave decay. The earlier Belle result which favored 1 ++ arose from the 
same sort of analysis applied to the 27r mass distribution in J/tpp. 

If this quantum number assignment is confirmed, it would remove the main motivation for the molecular 
interpretation of the A(3872), which assumes a loosely bound state of a D° and a D*° in S-wavc. It would also 
imply that the tetraquark hypothesis should be reexamined. Could it be, on the other hand, that the A (3872) 
is a standard 1 1 D2 charmonium? Beyond the known problems of its isospin violation in J/i/j p and J /%l> to [4], its 
larger than expected J / tpj and tj) 7 transitions [4-6] , and its high mass [7] , we would at least understand why 
the A(3872) is so narrow in the most prominent hadronic modes: they are close-to-threshold P-wave decays. 

We shall show (Sec. II) that there is another reason why the 1 1 D2 assignment is not so straightforward. Relying 
on an earlier study for D-wave charmonium fragmentation [8], we calculate the expected prompt production 
cross section at CDF and find that it is considerably lower than the observed A (3872) cross section; we also 
revisit the issue of the prompt production of the 1 ++ A (3872) in the molecular interpretation. As a 1 1 D2 
charmonium, the mass of the A(3872) is known to disagree with potential model predictions. We investigate 
here a different approach, a relativistic string model in the heavy quark limit (Sec. III). The model agrees 
remarkably well with the masses of established charmonia and bottomonia, but as in other approaches the 
A(3872) is very difficult to accommodate as a 1 1 D 2 state. We discuss possible decay modes that may help 
isolate the 1 X D2 assignment (Sec. IV). We also examine how to fit a 2 h A(3872) in a tetraquark model (Sec. 
V). The main drawback of this interpretation is the proliferation of predicted partner states, including many 
with masses close to the low- lying charmonia. We discuss the possibility (Sec. VI) that the A (3872) may be 
among the few which are narrow enough to be stable. 

The news about the A(3872) quantum numbers, if confirmed, will help discriminate among its possible 
interpretations. 



II. X(3872) PRODUCTION AT THE TEVATRON: MOLECULE AND CHARMONIUM 

In a recent paper [9] we proposed a method for estimating the prompt production cross section of A(3872) [10] 
at the Tevatron making the assumption that it is a loosely bound molecule of D° and D , with a binding energy 
as small as £q = —0.25 ± 0.40 MeV. The motivation for this study is that, after the Belle discovery [10], CDF 
and DO confirmed the A(3872) in proton- antiproton collisions [11, 12] and it seems at least counter-intuitive 
that such a loosely bound molecule could be produced promptly (i.e., not from B decay) in a high energy 
hadron collision environment. This was also one of the initial motivations to consider the possibility that the 
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A (3872) could be, instead of a molecule, a "pointlike" hadron resulting from the binding of a diquark and 
an antidiquark [13], following the interpretation proposed by Jaffe and Wilczek [14] of pentaquark baryons 
(antidiquark-antidiquark-quark). 

Indeed an analysis by CDF [15] allows one to distinguish the fraction of A(3872) produced promptly from 
the one originated from B-decays. The result of this analysis is 

a(pp -+ A(3872)) prompt x B(X -+ J/^ + ^) = (3.1 ± 0.7) nb. (1) 

Following [16], the experimental bounds on B(X — > J/ip-K + ir~) can be used to estimate a range for the prompt 
production cross section, 

a(pp -> A(3872)) prompt ~ 30 + 70 nb. (2) 

In Ref. [9] we estimated the maximum value of the theoretical prompt production cross section in the hypothesis 
that the A(3872) is a 1 ++ molecular state with a tiny binding energy So. Using standard MC tools, like HERWIG 
and PYTHIA, we computed the differential distribution of the prompt production cross section with respect to 
k, the modulus of the spatial component of the relative three momentum in the center of mass of the mesons 
constituting the molecular A (3872). This is a relevant variable for the system since we expect the two mesons 
to be almost collinear, otherwise the molecule would break down immediately (on a time scale r <~ 1/fc). In 
the MC simulation we have taken into account the kinematical cuts required for a comparison with CDF data. 
Integrating the differential distribution up to k < 50 MeV (an upper limit derived by the binding energy scale 
k w y/2m£ a ) we find that 

o(pp -> A(3872)) prompt < 0.1 nb, (3) 

about 300 times smaller than the measured value. This is a serious challenge to the molecular interpretation 
of the A (3872). 

On the other hand, in Ref. [16] it was argued that final state interactions (FSI) in the D°D*° system require 
two corrections to our previous calculation: (i) k should range up to A ~ 300 MeV; (ii) a correction factor 
to the cross section we compute should be included. This allows a spectacular reconcilement of the molecular 
picture with data, and Ref. [16] fully recovers the experimental prompt production cross section. 

In Ref. [17] we cast some doubts on the possibility that FSI can indeed play such a pivotal role. First, the 
Watson formulas [18] used in [16] are valid for S-wave scattering, whereas a relative three-momentum k of 
300 MeV indicates that higher partial waves should be taken into account. Most importantly, we have verified 
in our MC simulations that, as the relative momentum k in the center of mass of the molecule is taken to 
be up to 300 MeV, then other hadrons (on average more than two) have a relative momentum k < 100 MeV 
with respect to the D or the D* constituting the molecule. The "extra" hadrons are not only pions. On the 
other hand the Migdal- Watson theorem for FSI requires that only two particles in the final state participate 
in the strong interactions causing them to rescatter. The extra hadrons involved in the process interfere in an 
unknown way with the mesons assumed to rescatter into an A (3872). This is particularly true as one further 
enlarges the maximum value allowed for k, as required in [16]. 

In a recent paper [19] it is argued that the latter problem can be overcome as the interaction between the 
D mesons constituting the molecule is stronger than that between a £)W and one of the additional hadrons. 
Only a qualitative explanation of this argument is given. Moreover, the use of the MC approach to estimate 
the production of loosely bound hadron molecules proposed in [9] and then followed in [16] is criticized in [19]. 
We observe that: (i) The authors of [19] apply the method of [9] to deuteron production and, playing with the 
k parameter, they obtain disagreement with the empirical deuteron cross section by factors of 0(1), even if the 
MC generator is not tuned on baryon pair production data. We stress here that in the case of the A (3872) 
the MC cross section value is about 0.3% of the observed one, after having tuned the MC generator on D°D*~ 
data, (ii) The deuteron is a system qualitatively different from the DD* molecule, whose components cannot 
have spin interactions since the D is a spinless particle. In the case of the deuteron spin interactions play an 
important role in the determination of the binding: the spin S = deuterium is not formed. 

In our view the molecule picture is seriously challenged by the CDF cross section results. In any case if 
the 2 h quantum numbers are confirmed then the motivation for the molecular interpretation, which assumes 
S-wave binding, no longer applies. A 2 ^ state formed out of D°*D° would require a relative P-wave, and it is 
unlikely that tt exchange could bind such a state, given that even in S-wave it is not clear that the attraction is 
sufficiently strong [20]. Even if such a state exists, there remains the further problem that unless spin-dependent 
forces prevent the binding, one should expect partner states with h and the J PC -exotic 1 h state, for which 
there is no experimental evidence. A P-wave 2 ^ molecule would also imply the existence of a more deeply 
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bound S-wave 1 ++ molecule, which would be extremely narrow. Alternatively, to form a 2 h bound state in 
S-wave would require D 2 D or D\D* , which not only implies an immense binding energy of some 500 MeV, but 
one loses the appealing connection between the mass of the X(3872) and the D°*D° threshold. 

What about the prompt production rate if the X(3872) is a 1*D 2 standard charmonium? On general grounds 
one expects it to be small. The production cross section is proportional to a fragmentation function, which 
describes the probability that quarks and gluons hadronize into bound states. This function can be expressed 
as a perturbative expansion in the quark velocity v. For the production of a 1 1 D2 state the function begins at 
order 0(v 7 ), and so one expects a smaller cross section than for IP states, whose functions are 0(v 5 ), which 
are themselves suppressed with respect to IS states, 0(v 3 ). 

For the fragmentation functions we draw upon the result of Cho and Wise [8] , who calculated the produc- 
tion cross section of a 1 1 D 2 state at the Tevatron. They observed that despite the aforementioned kinematic 
suppression, large numerical prefactors in the amplitudes implied that such states could be produced in large 
enough measure to be observed in prompt production. These authors argue that color octet contributions are 
subleading with respect to color singlet in the heavy quark velocity expansion and thus are neglected. We 

use their result for the gluon fragmentation function D n „th)(z; /i) which describes the production of a *D 2 

9 ^r d 2 

quarkonia with helicity h, quark longitudinal momentum fraction z, and renormalisation scale /j,. We compute 
the production cross section, 

da 2 f 1 da 

-^(pp-^l 1 D 2 + AU) = V / dz-^-(pp ^g(p ± /z) + All; fi)xD (h) (z;n) (4) 
dp± t~o^° p± 

using recent gluon distribution functions. We find that the ratio between the MRSDO gluon distribution 
functions used in [8] and the most recent MSTW20008NNLO set amounts to about a factor of 0.7 in the 
most relevant Bjorken x region, x ~ y/s/s, which at Tevatron energies is x ~ V M 2 J 1960 2 ~ 0.02 , where 
M_l = ^JM 2 + p 2 _. Indeed we have M = 3872 MeV, p ± > 5 GcV, and \y\ < 0.6 to fulfil the kinematical cuts 
used in the CDF analysis. Here the factorization scale is set [i = M±. 

The integrated prompt cross section we find over the interval p± > 5 GeV is 

a(pp -> 1^2 + all) = 0.6 nb, (5) 

some 50 and 120 times smaller than the estimated experimental cross section in Eq. (2). It is difficult, therefore, 
to reconcile the observed production cross section of the A(3872) with the expectations for a 1 1 D 2 state. As 
we shall see, the mass of the X(3872) is also inconsistent with the 1 1 D 2 charmonium assignment. 



III. A RELATIVISTIC HADRON STRING MODEL IN THE HEAVY-QUARK LIMIT 

We turn our attention now to the question of the expected mass of a 1 1 D 2 state. Typically quark potential 
models predict a 1 X D 2 state lying some 50-100MeV lighter than the X(3872) mass [7]. In the following we will 
consider the mass of an 1 1 D 2 state in an evolution of the Chew-Frautschi string model (in the version described 
by Selem and Wilczek [21]) extended to accommodate heavy quarks. The string ends carry spin and electric 
charge, which allows us to include spin-orbit and magnetic moment interactions. 

As described in [21] (see also [22]), the energy of a relativistic string with tension T (=dE/dr in the rest 
frame of a segment dr) spinning with angular velocity oj is 

T r ri dv T r r2 dv 

£ = mm + to 272 + - / + - / (6) 

w Jo vl - v 2 w Jo \/l — v 2 

with 

y/l - (ojn) 2 

where ti is the distance of the particle i (i — 1,2) with mass rnj from the rotation axis. Analogously we can 
write the orbital angular momentum of the spinning string as 

2 2 t r ri v 2 dv t r r2 v 2 dv 

L = (7110)^71 + TO 2 wr 2 7 2 + — / + — / - (8) 

u 2 Jo VI - v 2 u 2 J Q VI - v 2 
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which derives from dL = ujdl = to (v 2 /uj 2 ) d£. Now, recall that the four-force is G^ = dp^ /ds, ds being the 
four-interval. Then G = 1/y/l — v 2 dp/dt = muj 2 j 2 r j, and j is in the centripetal direction. Thus we have 

T = mi"/ 2 uj 2 ri = m 2 "fiuj 2 r2 (9) 

Using Eq. (7) and (9), one can write -fi and as functions of T, m t and lj as 



li 



l ~\ 2 + 2 




(ii) 

Plugging (10) and (11) inside the expressions of £ and L one can obtain the energy and the orbital angular 
momentum of the system in terms of the mi, m 2 , T, and lj: 

£ = £{m 1 ,m 2 ,T,oj) (12) 
L = L(m\, m,2, T, oj) (13) 

In Appendix A we compute the expressions given above in the limit of equal and heavy masses mi = m 2 = M 
with T/Mlj « 1, or equivalently, TR « M. In this limit, 

E » 2M+^- (14) 



MTR 3 

L « ^— (15) 
with R = n + r 2 . In the heavy quark limit, therefore, the orbitally excited mesons have a fixed radius 

Miff (16) 

and there is a relation connecting the energy and the orbital angular momentum: 

£(L)*2M + 3(—") . (17) 

The string tension T is related to that of Regge phenomenology a by T = a/2n. This is in our picture the 
energy of a heavy quark-antiquark pair with L ^ orbital angular momentum. This is a characteristically 
different Regge trajectory from that of light quarks E <~ VtZ, which can also be obtained from the above 
parametric equations [21]. It has recently been demonstrated that the spectra of bottomonia does not fit with 
the standard Regge trajectory [23], and indeed we shall argue in this paper that the trajectory of Eq. (17) is 
the correct one to fit the data. 

We discussed this method in [24] and later found that the expression (17) had been previously mentioned 
in the literature [25], but not confronted directly with the meson spectrum. In Ref. [24] we applied it to the 
spectrum of orbitally excited tetraquark states in a diquark-antidiquark model. We note with interest that the 
same result, Eq. (17), arises in the Bohr-Sommerfcld quantization of the standard Cornell potential used in 
quark model calculations [26]. It shall come as no surprise, then, that the results presented here do not differ 
drastically from those of other models. 

We use Eq. (17) to predict the masses of the orbitally excited mesons with zero quark spin, yielding for the 
P- and D-wave spin singlets, 

M( 1 Pi) =2M + A, (18) 
M( 1 D 2 ) =2M + 2 2 / 3 A, (19) 

1 /3 

where A = 3 (T 2 /AM) ' . For the states with nonzero spin, we improve Eq. (17) by considering that at 
the ends of the string there is a heavy quark-antiquark pair carrying electric and color charge and magnetic 
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moment, which causes mass splittings due to spin-spin, spin-orbit, and tensor interactions. We thus write the 
mass formula for states with given total spin S, orbital angular momentum L, and total angular momentum J: 

M( 2S+1 Lj) = 2M + A L 2 ' 3 + M ss (2S q • Sq) + M LS (L ■ S) + M^(S 2 - 3 (5 • nf) (20) 

where here S q and Sq are the quark and antiquark spins coupled to S, and the unit vector n = R/R, with 
R the vector connecting the quark and the antiquark. The mass shifts M$s, Mls and M w for spin-spin, 
spin-orbit and tensor interactions are the analogues of the interaction potentials used in quark potential model 
calculations. The difference is that in this semiclassical approach, instead of taking matrix elements of these 
independent functions between meson wave functions we can evaluate the expectation values directly using the 
fixed value of the meson radius. We shall derive explicit expressions for these mass shifts later in the paper; 
first we investigate some general predictions of the model arising from Eq. (20). 

We assume that the spin-spin interaction splits only states with zero orbital angular momentum, for which 
the qq pair tends to be at much lower average distance. Assuming that M$s is common to both 1 So and 3 Si 
states, the masses of these states are 

MCSo) = 2M~^M SS (21) 

M( 3 Sr) = 2M+^M SS (22) 

Equations (18)-(19) and (21)-(22) yield a single linear relation among the masses 

1 - 2 2 / 3 

M( 1 D 2 ) = 2 2 / 3 M( 1 P 1 ) + (M (%) + 3M( 3 Si)) (23) 

and a parameter- independent prediction of the mass of the : D 2 state in terms of the masses of the h c , J/ip and 
Vc 

M( 1 D 2 ) ~ 3795 MeV. (24) 

The predicted mass is in reasonable agreement with other models, as we discuss later, but is considerably lower 
than the mass of the A (3872). In what follows we refine the model, adjusting the parameters to fit also the 
Xc states, but ultimately our predicted M( 1 D 2 ) changes very little. Thus, our conclusion is that the A(3872), 
whose mass exceeds our prediction by some 80 MeV, is difficult to reconcile with a 1 D 2 interpretation. 

In what follows we verify that the application of the same mass formula in the bottomonia sector is in 
remarkable agreement with data, supporting the validity of our prediction. It is not possible to apply the 
formula directly, since the bottomonia 1 Pi and X D 2 states have not been observed; instead we exploit the 
relation between the masses of these spin-singlet states and the spin-triplet members 3 Po,i,2 and 3 Di,2,3 of the 
same families. 

The spin-orbit and tensor terms split the states with spin-one and nonzero L, although we will assume, as 
is usual, that for these states spin-spin contact interactions are negligible. Applying Eq. (20) to the P-wave 
family yields the following expressions for the masses of the triplet states in terms of the mass of the singlet 
1 Pi and the shifts due to spin-orbit and tensor splittings: 

M( 3 P ) = M^Pi) -2M$ +2M<£> (25) 
M( 3 Pi) = MCP,) - M$ - M$ (26) 

M( 3 P 2 ) = MfV^ + M^ + l -M<f} (27) 

where we have labeled the mass shifts Mls and M MAi with a superscript to denote that they are dependent on 
the meson radius R, and hence the partial wave. Eliminating these terms yields a linear relation among the 
masses of the four members of the P-wave family, 

MCPx) = 1 (M( 3 P ) + 3M( 3 P0 + 5M( 3 P 2 )) . (28) 

The accuracy of this relation can be tested in the charmonia sector, where all four states have been observed. 
Using the Xco, Xci and Xc2 states as input, one predicts M{h c ) = 3525.30 ± 0.1 MeV, in remarkable agreement 
with the experimental mass M(h c ) = 3525. 67± 0.32 MeV [27]. Reference [28] uses recent high precision data on 
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the h c mass and finds that it is in even more striking agreement with this prediction. We note in passing that 
the good agreement also justifies the neglect of spin-spin interactions for nonzero L which, if present, would 
modify the relation (28). 

One naturally expects that the equivalent relation ought to hold in the bottomonia sector, and several authors 
have used it to predict the mass of the as-yet unobserved h b [29]. With the current masses of the \bo, Xbi and 
Xb2 we get 

M{h b ) = 9899.87 ± 0.59 MeV (29) 

Recent data [30] on the masses of the 3 Di i2 ,3 states of the upsilon sector allow us to apply the analogous 
approach, arriving at a new prediction for the mass of the as-yet unidentified 1 D 2 state. The mass formulas 

M( 3 D!) = M( 1 D 2 ) - 3M%> + MfP> (30) 

M( 3 D 2 ) = M( 1 D 2 ) - - MfJ (31) 

M( 3 D 3 ) = MCD 2 ) + 2M$ + |m#) (32) 



lead to the linear relation 



M( 1 D 2 ) = ^ (3M( 3 D!) + 5M( 3 D 2 ) + 7M( 3 D 3 )) (33) 



and we thus predict the mass of the unobserved *D 2 state, 

M(r) b2 ) = 10165.84 ± 1.8 MeV. (34) 

We are now in a position to test the validity of our mass formula (23) in the bottomonia sector. We use as 
input the experimental values for the r\ b and T masses, and for the h b , encouraged by the remarkable accuracy 
of the corresponding prediction in the charmonia sector, we use the center of gravity mass in Eq. (29). The 
mass relation (23) yields a predicted 1 D 2 mass of M(r) b 2) = 10168. 72+ 18 MeV, in striking agreement with the 
above value extracted from the experimental data, and therefore we can be confident in the reliability of the 
corresponding prediction for the 1 D 2 charmonia state. 

The remarkable accuracy of the model can be seen by plotting the Reggc trajectory of equation (17). For the 
center of gravity of the IS states £ (L = 0) = 2M we use Eqs. (21) and (22), 

£(L = 0) = J (3M(T) + Af (%)) = 9442.45+^° MeV (35) 

and for £(L = 1) and £ (L = 2) we use the values (29) and (34) extracted from the data. In the upper part of 
Fig. 1 we plot these three data points in the £-L 2 / 3 plane, and one sees immediately the accuracy with which 
the model fits the data. On the same plot (indicated by stars) we show our prediction for the masses of the 
higher L spin-singlet states, M( 1 F 3 ) = 10391 MeV and M( 1 G 4 ) = 10591 MeV, the identification of which would 
be a good test of our model. Both of these states should be narrow since they lie below the threshold for B*B, 
the lightest open flavor pair to which they can decay. 

In the lower part of the same figure we plot the corresponding Rcgge trajectory for the cc system. Here we 
have only two data points: £(L = 0), which we determine from the masses of the J/ip and r\ c as above, and 
£{L = lj = M(h c ). We indicate by stars our predictions for M( 1 D 2 ) = 3795 MeV, M( 1 F 3 ) = 4020 MeV and 
M( 1 G 4 ) = 4221 MeV, although the masses of the latter pair may not be as reliable as in the bottomonia case 
because of the coupling to the open D*D threshold. 

A striking feature of Fig. 1 is that the slopes of the trajectories for bb and cc are almost identical. Since the 
slope is given by 3(T 2 /4M) 1 / 3 , this implies that the effective string tension scales with the square root of the 
quark mass, in contrast to many potential models in which a common string tension is used for different quark 
masses; we determine the best fit values T b = 0.257 GeV 2 and T c = 0.148 GeV 2 . From the intercept of the 
trajectories we determine the quark masses M b = 4721 MeV and M c = 1533 MeV. 

We return now to the derivation of the mass shifts M$s, Mls, and M M/J in terms of the parameters of the 
string model. In a potential model approach the interaction potential for the spin-spin term is 

V SS (R) = k c5 5 3 (R). (36) 

The corresponding mass shift is thus proportional to the square of the wave function at the origin, assuming 
that the magnetic coupling and the wave function at the origin do not depend on the total spin of the pair 

M ss = (V SS (R)) = i^cc |^s(0)| 2 • (37) 
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FIG. 1: Regge trajectories for bottomonia (upper left) and charmonia (lower right) states. The closed circles indicate 
experimental masses (or those of the center of gravity where appropriate) for a given L; stars indicate predictions for 
unobserved states. 



In our model we treat k' = k 1^(0)1 2 as a free parameter. The mass shift Mls is due to the coupling of the 
magnetic moment of one quark to the magnetic field created by the moving charge of the other quark; this is 
the essence of the spin-orbit coupling, 

where A is a constant to be determined and £ is the energy of the string as in Eq. (6) with R = r% + r 2 , an 
approximately linear function as shown in Eq. (14). The tensor mass shift arises from the interaction 

between the magnetic moment of one of the two quarks with the static magnetic field generated by the other 
one, 

M pp = B± (39) 

where the constant B is again to be determined. 

Following standard electromagnetism, one can parametrize A and B in terms of the gyromagnetic factor and 
of the charge of the heavy quarks bound to form the meson. The interaction between the magnetic dipole and 
the magnetic field generated from a moving charge carries a factor ge/2m 2 , whereas the Thomas precession 
gives — e/2m 2 . As for the tensor coupling B we expect it to be related to the product of the Bohr magnetons 
of the quark and the antiquark. Thus, one is led to write 

(ff Q ~ 1 )' Q _ / -v...' ,j 



As we are dealing with a constituent quark model we do not expect the gQ to be comparable to those of a 
pointlike structureless particle. Thus, we leave gQ as a free parameter and, as expected, we find best fit values 
larger than 2. We use e c = 2e/3 and e& = e/3. 

Wc test the validity of our model on the charmonia and bottomonia spectra with L = 1, 2. In the heavy quark 
limit we expect the model to work better at fitting the bottomonia states, and we fit the rji,, T, Xbo, Xbi,Xb2 [27], 
and T(lDi), Y(1D 2 ), T(1_D 3 ) [30]. The splitting between rjb and T is not a feature of the string model per se, 
since it is controlled by the parameter n' h which is not correlated with any other masses. The masses of the 
remaining states are controlled by two parameters (Tf, and gb), and we find the best fit values 




0.258 GeV 2 



(41) 
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FIG. 2: Results of the fit to the L = 1 bottomonia FIG. 3: Results of the fit to the L = 2 bottomo- 

with the parameters in Eq. (41). nium states with the parameters in Eq. (41). 



The string tension is close to the value used to fit Regge trajectories of light mesons (T w 0.175 GeV 2 ) and differs 
only very slightly from the value obtained earlier by fitting the spin-averaged masses. The resulting spectrum is 
summarized in Table I, where we also compare our results with the predictions obtained in the potential models 
of Refs. [31] and [32]. The agreement is remarkable. The splitting between the center of gravities of the S-, P-, 
and D-wave sectors has already been confronted with data; the new feature is that with the same parameters 
we are also able to describe the splitting among the 3 Po,i,2 and 3 Di2,3 states. A comparison of the theoretical 
and experimental spectra is presented in Figs. 2 and 3. 



M exp (MeV) 


a(M exp ) (MeV) 


M th (MeV) 


M th (MeV) [31] 


M th (MeV) [32] 


Vb 


9391 


3.2 


9388 


9366 


9400 


T(15) 


9460 


0.26 


9460 


9460 


9460 


h b 




9898 


9924 


9880 


Xbo 


9859 


0.5 


9861 


9888 


9850 


Xbl 


9893 


0.4 


9900 


9913 


9880 


Xb2 


9912 


0.4 


9905 


9939 


9900 


T(1 1 D 2 ) 




10166 


10166 


10150 


T(l 3 Di) 


10152 


1 


10147 


10153 


10140 


T(1 3 D 2 ) 


10164 


0.9 


10163 


10163 


10150 


T(1 3 D 3 ) 


10173 


1 


10177 


10174 


10160 



TABLE I: Theoretical values of the masses obtained by the fit with the parameters in Eq. (41) compared to the experi- 
mental values and to other theoretical determinations [31, 32]. 

The same work can be done for charmonia. In this case a priori we do not expect such a remarkable agreement 
with data. Still, as we will see, the agreement with the L — 1 states is rather good. We include in our fit the 
following charmonia [27]: r) c , J/ip, h c , XcO>Xci>Xc2> and ^(3770), which we identify with the l 3 Di state. The 
best fit is obtained with the following values for the parameters: 

T c = 0.147 GeV 2 
9c = 5.7 

Here the agreement of the value obtained for the string tension with the one used for light mesons is better 
than in the bb case, and again it differs very little from the value obtained in the spin-averaged case. It is worth 
saying that performing the same fit procedure numerically with the full expression for the energy, i.e., without 
any expansion, and orbital momentum [see Eqs. (6) and (8)], the best fit parameters are essentially the same 
as the ones quoted above, both for the charmonium and the bottomonium sectors. 

The experimental masses [27] are compared to the results of the fit in Table II, where we again compare with 
the predictions of Refs. [31, 32], and in Figs. 4 and 5. The agreement is very good, although unlike Refs. [31, 32] 
our model fails to reproduce the mass of the ^(3770), which is around 45 MeV above our predicted value. The 
effect of the mixing with the nearby i//(3686) could be advocated as a possible resolution of this mismatch. 
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FIG. 4: Results of the fit to the L = I charmo- 
nium states with the parameters in Eq. (42) . The 
agreement between data and the results of our 
string model are remarkable in the L — 1 sector. 



FIG. 5: The same as in Fig. 4 but for the L = 2 
states. We tentatively identify the X(3872) with 
the l 1 r>2 charmonium state and the 0(3770) with 
the l 3 Di. 



However, since the mixing angle between the two states is only <j> = (12 ± 2)° [33, 34], it could account for an 
upward shift of at most ~ 5 MeV with respect to the fitted value, giving the -0(3770) at 3735 MeV. 



M exp (MeV) 


a(M exp ) (MeV) 


M th (MeV) 


M th (MeV) [31] 


Mth (MeV) [32] 


Vc 


2980 


1 


2981 


2965 


2970 


J/0 


3096 


0.1 


3097 


3095 


3100 


h c 


3526 


0.32 


3525 


3525 


3520 


XcO 


3415 


0.31 


3427 


3415 


3440 


Xcl 


3510 


0.07 


3518 


3508 


3510 


Xc2 


3556 


0.09 


3548 


3555 


3550 


V(3872)[1 1 D 2 ] 


3872 


0.24 


3793 


3810 


3840 


0(3770) [l 3 Di] 


3772 


0.35 


3727 


3762 


3820 


1 3 D 2 




3779 


3797 


3840 






3831 


3840 


3850 



TABLE II: Theoretical values of the masses obtained by the fit with the parameters in Eq. (42) compared to the 
experimental values and to other theoretical determinations [31, 32]. 

Our predicted value of 3793 MeV for the mass of the 1 1 D2 state is only slightly different from that which we 
obtained using the formula (24), and our conclusion remains the same: the mass of the AT (3872) is difficult to 
reconcile with the charmonium interpretation. Indeed it is evident in Figs. 4 and 5 the most difficult identification 
to make is that of the AT(3872) as a 1 X D2 state. As can be seen in Table II, the same is also true of the potential 
models [31] and [32]. Reference [7] compiles the mass predictions for a variety of different potential models, and 
five out of the six cases (including the two which we have quoted in this paper) predict the mass of the 1 1 D2 is 
some 50-100 MeV lighter than that of the X(3872). The exception is that of Fulcher [35] who predicted a 1 1 D2 
with a mass of none other than 3872 MeV; notwithstanding the remarkable agreement of that model with the 
AT (3872) mass, we find rather more compelling the broad agreement of our model with the predictions of the 
majority of other approaches. We note in passing, however, that lattice QCD predicts a somewhat higher 2 h 
mass of 3907 ± 32 MeV [36], consistent with the A" (3872). 



IV. DECAYS OF : D 2 CHARMONIUM 



We have discussed two difficulties with the charmonium interpretation of the ^(3872): its prompt production 
is too small, while its mass is too large. Another difficulty is the isospin violation in its decays. Indeed 
the J/tpp and J/ipuJ modes are observed with the same strength: B(X — » J/ipu)/B(X — > J/ipTr + ir~) = 
1.0 ± 0.4 (stat.) ± 0.3 (syst.) [4], in obvious contradiction to expectations for a standard cc state which is a 
pure isoscalar. However the isospin violation is not so severe if one considers the different phase space volumes 
available to J/^tt + tt~ and J/07r + 7r~7r° final states. Because of the different decay widths of p and w one finds 
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the ratio between the I — 1 and 1 = amplitudes Ai=i/Ai = o ~ 0.5, consistent with the experimental value. 
A further possibility is that the J/07r + 7r _ and J / tpTT + n a modes are fed by rescattering from intermediate 
D*D states, and isospin is broken by the mass difference between the neutral and charged states [37]. 

Historically, the literature on the 1 1 T>2 state has assumed that it lies below the D*D threshold in which case 
it would decay dominantly by radiative transitions, hadronic transitions with pion emission, and annihilation 
into gluons. Among the radiative transitions, the dominant mode is expected to be [8] 

A(3872) -> h cl -> J/V>7T°7, (43) 

with a branching ratio of 0.004. Using the A (3872) mass, Ref. [7] predicts a h c j width of 0.460 MeV, somewhat 
larger than an earlier prediction of 0.278 MeV which assumed a lower mass [38]. Reference [7] also predicts 
a considerable -0(3770)7 mode with a partial width of 0.045 MeV. By contrast, the observed J/V>7 and ip'j 
are rather at odds with expectations. Because of the orthogonality of the spin and spatial wave functions 
of the heavy quark pair, these transitions are expected to be small, although the V'7 would be enhanced 
due to its ID component, as it is presumably the orthogonal partner of the dominantly ID ■0(3770). In a 
recent paper [6] the 2 h hypothesis has been tested on the radiative decays of A (3872). Using potential non- 
relativistic QCD (coupled to electromagnetism to describe single photon transitions), they found upper bounds 
on i3(l 1 Z?2 — > J )j) which are 1 order of magnitude smaller than the experimental lower bounds obtained 
from BaBar measurements. The estimate for the 1 1 I?2 — > "0 7 takes into account the S — D wave mixing in the 
i\) wave function, which turns out to be not sufficient to accommodate theoretical predictions with experimental 
data. The approach used in [6] , first derived in [39] , relies on the use of charmonium potential models in order 
to compute overlap integrals of radial wave functions. The stability of the results is checked using different 
potentials. Furthermore in recent news Belle does not confirm the A (3872) — > tp'j decay [40], in contrast to the 
BaBar result with an estimated branching ratio B(X — >• ip 7) ~ 6% [41]. 

As for the hadronic transitions Refs. [7, 38] indicate that the dominant hadronic decay mode with two pions 
in the final state is 1 1 D2 —> rj c Tnr, and the latter authors predict a partial width of 0.21 ± 0.11 MeV. As 
described in [42], this result is obtained using a multipole expansion of the color gauge field. In atomic physics 
the analogous process is the emission of atomic electrons or e + e~ pairs by nuclei. This expansion leads to some 
selection rules which at leading order imply that quark spin is conserved. Thus the operator responsible for the 
transition connects quarkonium states with the same spin S, namely 1 1 D 2 and l 1 So, the r\ c . The observation 
or the absence of this decay would be crucial and thus we encourage the experimental search in this direction. 

Reference [43] considered the D°*D° decays of a heavier 1 1 D 2 and predict a partial width of 0.03 MeV for a 
state with mass 3872 MeV, increasing to 1.7 MeV if it were somewhat higher at 3880 MeV; experimentally, the 
A (3872) is known to lie much closer to 3872 MeV and these partial widths are consistent with the experimental 
width. 



V. THE TETRAQUARK INTERPRETATION 

In this section we consider the tetraquark option which, unlike the molecule, may still be a viable interpretation 
of the A (3872). One advantage of the tetraquark interpretation is that it admits a natural solution to the isospin 
violation issue, due to suitable mixing between the diquark flavor states [13]. To form a diquark-antidiquark 
state with 2 h quantum numbers we need a unit of orbital angular momentum between the two diquarks (the 
diquarks have positive parity), and at least one out of the two diquarks must have spin 1. The existence of 
such a state therefore implies the existence of a rich spectrum of partner states, formed of combinations of spin 
1 (axial vector) and spin (scalar) diquarks coupled to the L = 1 orbital angular momentum to give various 
total angular momenta J; in addition, one expects a series of lighter L = states with positive parity. The spin 
wave functions and corresponding J PC of L = and L = 1 tctraquarks is discussed in Ref. [44]; we extend 
their results to states with J = 2 and 3 and summarize these in Table III. We shall assume, as is usual, that 
only color (anti)triplet diquarks form stable configurations [45]. 

In general a tetraquark with a given J PC can be formed from one of several different spin wave functions, and 
we thus expect that the physical states are admixtures of the various possible states of (S cq x S S q)s- Fot a 2~+ 
tetraquark there are two possibilities: an orbital excitation of a vector diquark pair coupled antisymmetrically 
in spin or a scalar and vector diquark likewise coupled antisymmetrically. One needs also to consider the various 
flavor combinations; there should exist, for each angular momentum configuration, two neutral and two charged 
states formed out of [cu] and [cd] building blocks and the corresponding antidiquarks (and a further four strange 
states if one considers strange diquarks) . This presents something of a problem from a phenomenological point 
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TABLE III: The spectrum of L = and L = 1 tetraquark states classified according to their spin wave functions (with 
Sea, Scs, and S the diquark, antidiquark, and total quark spin, respectively) and J PC quantum numbers. 

of view; if the V(3872) is a 2~ + tetraquark state then, a priori at least, it should be accompanied by at least 
seven partner states, three neutral and four charged. 

It is possible that the data reflects the existence of two neutral states, the lighter state decaying to J/ipp and 
J/ipuj and the heavier to D°*D°, although the experimental interpretation is not clear [46, 47]. In Ref. [48] 
this possibility was discussed in the context of the X(3872) as a 1 ++ tetraquark, where a total of two neutral 
states and two charged states are expected. Even if there are two states in the data, in the 2 h interpretation 
there would still be two missing neutral states. Moreover, we expect four charged states, for which there is 
apparently no experimental evidence [49] . The problem may not be so severe if one takes account of the mass 
difference between scalar and vector diquarks, which in the relativistic approach of Ref. [50] is around 60MeV. 
In this case the physical X(3872) could be identified with the state dominated by the scalar- vector combination 
in Table III, and the vector- vector combination would be somewhat heavier. If the P-wave coupling to J/^lo, 
J/^p, and D°*D° blows up rapidly with momentum, then this heavier state may be too broad to be observed; 
in this case one still needs two neutral and two charged states, instead of four and four. 

We also expect the X(3872) to be accompanied by a staggering array of partners with different J PC : one 
each of 0~~ and 3~~, two each of 2 , '", 1 h and 2 and four For each J PC we again expect two 

neutral and two charged members, and possibly also strange states. If the S-wave states also exist one expects, 
also multiplied fourfold in flavor, two each of ++ and l + ~, and one each of 1 ++ and 2 ++ . In the absence of 
experimental evidence for these partner states a tetraquark explanation for the AT (3872) is only tenable if one 
can explain why it, among all the possible configurations, is unique. We shall propose one possibility which 
is that most of the states are so broad as to be effectively unobservable, while the AT (3872), due the scarcity 
of allowed hadronic decays and limited phase space, may be among the few states which are narrow and thus 
observable. 

Spectrum. Before addressing such questions we consider the masses we expect for these tetraquark states. 
For the L = 1 states there will be presumably be some spin-orbit splittings, although we do not calculate 
these here. We shall assume, instead, that the states lie around the mass of the V(3872), somewhere in the 
3800 — 3900 MeV mass region. Moreover, we assume that the scalar and vector diquarks have equal masses. We 
fit the diquark mass with the relativistic string model (17) in such a way that at 3872 MeV we find the L = 1 
state. This implies 

m [cg] = 1716 MeV, (44) 

where we use for the string tension T the value obtained from the fit to the cc spectrum. In Ref. [13] we 
discussed the interpretation of the A (3872) as an S-wavc 1 ++ state, and correspondingly the diquark mass 
mu] = 1933 MeV was larger than what we obtain here. 

The entire spectrum of L = states discussed in [13] is shifted downwards in mass once the new input value 
is the 2 h V(3872): the average mass will be around 2mr c ,i ~ 3430 MeV. Thus a striking prediction of the 
tetraquark model is that there is a set of positive parity states in the same mass region as the Xco.1,2 and h c . 
Their masses can be expressed in terms of the diquark mass and of the magnetic couplings among the four 
quarks which constitute the state, since at L = one can assume that only spin-spin interactions are relevant 
[13]. The full spectrum of tetraquark states is reported in Fig. 6. Among the L = states the ++ tetraquark 
is very close in mass to the standard charmonium state Xco, while the others are some 40 MeV lighter than the 
charmonia with the same J PC . 

As for the L — 1 part of the spectrum, a closer look is required for the 1 states. In a previous paper [24] we 
proposed a tetraquark interpretation for the vector resonances V(4350) and Yb (4660). Motivated by a reanalysis 
of the Belle data on the Yb (4660) which shows a preference for the baryon-antibaryon decay mode (the most 
striking tetraquark signature), these states were identified with L = 1 tetraquark states with zero and one unit 
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FIG. 6: Dashed lines represents the full spectrum of the tetraquark [eg] [eg] states with L — and L = 1, whereas the 
solid lines are the standard charmonium levels with L = 1. 



of radial excitation, respectively. The masses of the states were obtained using Eq. (17) with M — 1933 MeV, 
the value obtained for the diquark mass in [13] in the hypothesis that the X(3872) were a 1 ++ state. Reference 
[50] considered a similar scheme in a relativistic diquark-antidiquark model and likewise obtained a spectrum 
of L — 1 states in the 4200-4300 MeV mass region. 

The new value of the diquark mass (44) pushes the masses of the L = 1 states down to the 3800 — 3900 MeV 
region. In this scheme the Y(4350) and the Yb (4660) could thus be the two successive radial excitations of the 
l -- state at about 3870 MeV . For a rough estimate of the mass splitting between different radial excitations, 
we borrow the results obtained for the P-wave charmonia [51], which indicate a mass gap from the ground to 
first radial excitation of around 440 MeV, and from the first to second radial excitation of around 370 MeV. 
Using this approach, as in [24], we predict 

M(n r =1,L = 1) = 4320 MeV 

M(n r = 2, L = 1) = 4690 MeV ^ ' 

The advantage of this assignment is that it explains the suppression of the Y(4350) — > J j^m^ix^ and 
Yb (4660) — > J/ijm + TT~~ decays with respect to the corresponding tfi decay modes. The J/ip modes should 
be enhanced by phase space, but may be suppressed due to different radial quantum numbers between the 
charm quarks in the initial and final states. In [24] it is shown that the measured spectra indicate that the 7T7t 
pair comes from phase space (or, equivalently a a) apart from the case of the Yb (4660) where there seems to 
be a 30% component due to /o(980). As it was done in [24], we attempt an explanation of the observed Y's 
decay pattern describing the S-wave transition (ip(lS, 25)a|Y) as 

(iP(r,,q)a(k)\Y(e,p))=ge-r) (46) 

where a = cr, fo- Interpreting the Y's as charmoniumlike bound states made up by a diquark and an antidiquark, 
g can be estimated from the radial overlap integral, 

5 cx J d 3 r R^R^ir) (47) 



1 Observe that using the tension obtained from the fit to the charmonium states, the mass of L = 1 and L = 3 tetraquarks would 
be 3872 MeV and 4353 MeV respectively. We choose not to consider the L = 3 assignment. 
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where R (Y) (r) is the radial wave function for the P-wave [eg] [eg] states and i?^(r) refers to the standard 
cc states. Since our string model provides us with an expression of the interaction energy as a function of 
the distance between the quark and the antiquark in the case of a charmonium state, and the diquark and 
the antidiquark in the case of a tetraquark, we solve numerically the radial part of the Schrodinger equation 
using for the potential equation (A3) with the parameters T and M obtained from the charmonia fits. We are 
interested in estimating the following ratio of decay widths: 

r a (F) = T(Y -+ ^(lS)(7T7r) a )/r(y -> ^'(7r7r) a ) (48) 

The square of the matrix element weights the three-body phase space where the two pions are the decay 
products of an intermediate scalar resonance. The Breit-Wigner ansatz we use is a rough approximation for 
the description of a very broad a meson whereas is rather suitable for the /o(980). We find for the Y(4350) 
the ratio T a = 0.3, to be compared to the limit extracted in [24]: £(Y(4350) -> J/^tt + tt-)/B(Y(4350) -> 
ip'ir+ir-) < 3.4 x 1CT 3 @ 90 %C.L. For the Y B (4660) the corresponding ratios due to both a and / (980) 
modes are T CT = 0.6 and T fo = 1.9, to be compared with #(Y B (4660) -> J/ipn + 'K-)/B(Y B (A660) -> ip'ir+ir-) < 
0.46 @ 90 %C.L.[24]. Although it cannot reproduce entirely the experimental suppression of the J/ipir + n~ 
modes, the ansatz contained in Eq. (47) is able to account for some enhancement of tttt despite the phase 
space suppression. 



VI. DECAYS OF TETRAQUARK STATES 

We give now some general considerations about the expected decay patterns of these tetraquark states. Each 
tetraquark wave function contains an admixture of color singlet meson pairs, and one therefore expects that 
hadronic decays will be dominated by a "fall-apart" mechanism in which the tetraquark dissociates into meson 
pairs, either charmed mesons or a charmonium plus light meson(s), and also radiative decays in which a light 
quark pair annihilates into a photon. The A (3872) decays in all three ways. 

We show in Table IV the J PC, -allowed hadronic and radiative decays of L = and L = 1 tctraquarks and 
the corresponding partial waves. We also identify selection rules which arise due to the spin part of the decay 
amplitude, the details of which are derived in Appendix C. Entries without parentheses indicate that the 
fall-apart decay can proceed by a spin-conserving process, which is to say that the diquark-antidiquark spin 
wave function recouples directly into the spin wave functions of the final state mesons. Those with parentheses 
require the nonconservation of spin, either of the light quark or heavy-light pair (round brackets), or the heavy 
quark pair (square brackets); we expect the latter to be a stronger selection rule. Such arguments have been 
used before in the context of the X(3872) as a 1 ++ tetraquark [13] or molecule [52], where the wave function in 
either case necessarily has the cc in spin-one, consistent with the dominance of the J/-07r + 7r" and J / ipir + ir~ ir° 
modes. 

The L—0 states. The J + states present something of a problem from a phenomenological point of view, 
insofar as they are as light or even lighter than the corresponding h c and \c states, sharing the same quantum 
numbers and hence decay modes, and yet there is apparently no experimental evidence for their existence. Recall 
that among the neutral states we would expect four additional ++ and 1 H , and two additional 1 ++ and 2 ++ 
states. Referring to Table IV, we note that the two most numerous of the expected additional states, ++ and 
l" 1 , can decay in S-wave to r) c Tr and J/i^tt respectively, and one could argue that they simply "fall-apart" 
broadly in such a way as to be effectively unobservable. The remaining states L — cannot be dismissed in this 
way. The 1 ++ states, because of their unnatural parity, cannot decay into t] c tt and we thus expect them to be 
comparatively narrow. The corresponding Xd state has a width of less than 1 MeV; the tetraquark analogue 
may be broader on account of fall-apart decay into r] c irir via the low mass tail of the broad a, however this is 
forbidden in the limit of heavy quark spin conservation, since the 1 ++ tetraquark necessarily has the cc in spin 
1, as discussed above. The 2 ++ states could decay to r] c Tr in D-wave, but this is also forbidden by heavy quark 
spin conservation. In the tetraquark picture we thus expect light, narrow 1 ++ and 2 ++ states in the Xc mass 
region decaying into r/ c 7r and J/V>7- 

The L—l states. Because of their higher mass the J~ states have many more available decay modes and 
considerably more phase space, and it is plausible that most of them decay broadly so as to be effectively 
unobservable. The A (3872) with 2 ^ quantum numbers may be a unique exception. Because of its unnatural 
parity it cannot decay to i] c tt 7 r] c r], or DD; its observed decays into J/ipp, J/ipuj and D°*D° are all P-wave 
with very little phase space, naturally implying a small width in accordance with the experimental data, T = 
3.0^'y MeV. Moreover, we note that the J/ip modes are suppressed by the conservation of quark spin, although 
because only the light quark spin needs to flip this selection rule may be badly violated. Indeed the very 
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TABLE IV: Hadronic and radiative decays of tetraquark states in given partial waves S, P, D, F. Partial waves in 
parentheses indicate that the decay can take place only by spin-flip, either of the light or heavy-light quark pair (round 
brackets) or the heavy quark pair (square brackets). In applying the rule to decays involving the a, we assume that they 
proceed through an intermediate qq with 3 Po quantum numbers which then feeds the physical a. 



observation of these modes confirms that the rule is broken, but without further model-dependent assumptions 
it is difficult to ascertain to what extent the narrowness of the X(3872) is due to the limited phase space and to 
what extent it is due to this weak selection rule. We note also that the observed J/ tp^f and tp 7 decays also imply 
the nonconscrvation of light quark spin: the qq pair can only annihilate into a photon if it has spin one, and the 
wave function of the X(3872) does not contain such a component (see Appendix C). In a diquark-antidiquark 
model the P-wave tetraquark wave function contains a cc pair at large distance due to the P-wave, and thus 
it may be anticipated that Xc^ modes, particularly Xc2^, and h c -f, should be large. We thus urge a search for 
these challenging modes which, if observed to be prominent, would support the tetraquark hypothesis, although 
as noted earlier h c j should also be large in the 1 1 D2 case. 

The pattern of allowed decays for the 0~ + states is very similar to those of the 2 h states, largely because of 
their shared unnatural parity. The r\ c a decay goes in S- rather than D-wave, and in a model approach [13] one 
expects tetraquark states to couple strongly to final states which themselves have four-quark content, as the a 
is generally accepted to have. This coupling may be so strong as to make the width of the h much larger 
than that of the 2 which may therefore make it more difficult to identify. If such a mechanism is not in place 
then the ^ states should exist and will be as narrow as the AT (3872). Indeed, they may be even narrower: 
due to spin-orbit splittings one expects the ^ to be lighter than the 2 h [44], which would close the D°*D° 
mode and reduce the effective phase space for J/tpp and J/ipw. Such a state should, like the A'(3872), decay 
radiatively into J/ip-y, and there is apparently no signal in the data [41]. 

Among the remaining states we expect the unnatural parity 0~~ and 2~~ to be the most narrow, although 
they will presumably be broader than the corresponding ^ and 2 ^ states due to the much greater phase 
space available to J/tpir and J/xp-q compared to J/^pp and J/tpu; if the decay amplitudes blow up rapidly as 
phase space opens up, these may be too unstable to be identified. If their mass is sufficient the D*D mode could 
be prominent, while r\ c p and r] c uj are forbidden by the conservation of heavy quark spin. We expect dominant 
radiative decays to XciJ and, for the 2~, Xc27- 

The remaining 1~~, 3~~, and J PC -cxotic 1 _+ states have many decay modes available and with ample phase 
space, so we expect that they should be the least stable of all the possible L — 1 tetraquark configurations. 
This may be helpful from a phenomenological point of view. The spectra of 1~~ states in the 3800-3900 McV 
mass region has been very well studied and there is apparently no evidence of an overpopulation of states, 
whereas in the tetraquark picture one expects a further eight neutral states alone. The phenomenology is thus 
only consistent with the assumption that the 1 states, due to the large number of decay modes available and 
ample phase space, do not exist as stable resonances. It may be, however, that their radial excitations are more 
stable on account of the spatial separation of the cc pair, as discussed earlier in the context of the y(4350) 
and Yb (4660). If it is indeed the case that the radial ground states of the L = 1 tetraquark states with 1 
quantum numbers do not exist as stable resonances, then one can infer, judging by the overall similarity of their 
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decay patterns to those of the 1 states, that the 3 and 1 + states are likewise unstable and probably do 
not exist. 



VII. CONCLUSIONS 



The BaBar result which favors 2 h quantum numbers for the A(3872) implies a serious revision of theoretical 
interpretations is required. If it is confirmed, the molecular interpretation appears to be untenable, but as we 
have shown the 1 1 D 2 and tetraquark interpretations also face challenges. 

If the A (3872) is a 1 1 D2 state, we would expect a much smaller prompt production cross section than that 
which is observed at the Tevatron. Moreover, we have confirmed the result obtained in other models that the 
mass is much heavier than would be expected. Our model, an extension of one previously discussed in [24], 
predicts with remarkable accuracy the masses of most established charmonia and bottomonia mesons, while for 
the X(3872) there is a considerable discrepancy. As for the decay pattern we suggest that h c j and r] c irir ought 
to be prominent, and that the verification or otherwise of 1//7 may help clarify the situation. 

We reexamined the tetraquark spectrum in the hypothesis that the A(3872) is a P-wave tetraquark. The 
corresponding diquark mass is lower than that obtained if A(3872) were a 1 ++ state, implying a series of S- 
wave states with masses comparable to P-wave charmonia. This may be difficult to reconcile with data, and, 
in particular, we expect narrow 1 ++ and 2 ++ states. On the other hand, among the P-wave tetraquarks in the 
3800 — 3900 MeV region it is the A(3872) with 2~ + quantum numbers which ought to be the most stable; many 
of its partner states, for which there is no experimental evidence, are probably too broad to exist as genuine 
resonances. On the other hand, we must emphasize that the narrowness of the A(3872) is not a consequence of 
the assumed four-quark structure: any 2 h state at 3872 MeV will be narrow because of the remarkable fact 
that its allowed decays are all P-wave with almost no phase space. An important consequence of the tetraquark 
interpretation is the reassignment of the F(4350) and Yb (4660) with respect to Ref. [24]; in this picture they 
could be understood as the first and second radial excitations of a P-wave [cq] [cq] bound state, explaining why 
the tjj 7T7T decay mode dominates over the J/ipnir one despite of phase space suppression. 



Appendix A: Equal and Infinite heavy mass limit 



The energy and the orbital angular momentum of a relativistic string spinning with an angular velocity ui 
with two masses m\ and m 2 attached at each end can be written as a function of T, mi, m 2 and oj. In the 
equal and infinite mass limit (mi = m 2 = M), we can expand Eq. (12) up to the second order in the parameter 

e = T/Moj, 



2M 



.2T 



e 2 M + C(e 3 ) - 2M 



3T Z 
Mu 2 



(Al) 



IT 2 - 
Muj 3 



In this limit we can compute the dependence of £ on the distance between the masses, using Eq. (9) to 
eliminate co in favor of R = 2r = r\j2 = r 2 /2: 



One finds 



uj 2 = (1 - (Lur) 2 )^- ^> lu = J - - 9rTi (A2) 
v v ' ' Mr V Mr + r 2 T 



£ =2M+^f + ^f 

3/2 (A3) 



T 2 / 2MR+R 2 T \ 
4M V T / 

which reduce to Eqs. (14) and (15) for M » TR. From the first of these two relations it is evident that this 
relativistic string model well describes the confinement, the energy being a growing function of the distance 
between the quark and the antiquark inside the meson. The second of the relations in (A3) allows to write the 
distance R between the quarks in terms of the orbital angular momentum, the mass, and the tension: 

M ^fM 2 T 2 /z + (4LMT 2 ) 2 / 3 

R -~Y + T^T 3 ( ' 
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which reduces to Eq. (16) for M » TR. 



Appendix B: Relativistic Corrections 

The spin-orbit term contains the derivative of the energy E with respect to the distance R. Using the first 
relation in Eq. (A3) gives 

dE 3T(M + RT) 

dR 2M (B1) 

The spin orbit term can thus be rewritten in terms of T, M and L plugging the expressions found for R in (A4) 
inside Eq. (Bl). For the tensor term the magnetic field generated from a magnetic moment Mi is 

3n (pi -n) - mi fno \ 

Hi = ~ 3 (B2) 

The interaction energy is 

M w = . Hl -mi H 2 = _ 3(/x 2 -n) (j^-j^ + (1 ^ 2) ^ 

where n is the unit vector in the direction between the two heavy quarks. With 

m - (B4) 

and 

S 2 - 3 (5 • nf = 2 [Si ■ 5 2 - 3 (5i • n) (5 2 • n)] (B5) 

one finds 

[5 2 -3(5-n) 2 j (B6) 
Again, using the expression (A4) we can write the tensor interaction as a function of T, M and L. 



Appendix C: Spin recoupling matrix elements 

The diquark-antidiquark spin wavefunctions in the (S cq x S B q)s basis can be rewritten in the bases (S C q x S q5 )s 
and (S cS x S q q)s of the open- and closed-flavor final states, respectively. The matrix elements ((S cS x S q q)s\ (S cq x 
Scq)s) f° r closed flavor decay are 

|(0x0)o) |(lxl)o) 
((0x0)o| 1/2 V3/2 (CI) 
((1 x 1) | V3/2 -1/2 



(C2) 



1(1x1)!) 0/2(1(1x0)0 + 1(0x1)0) 0/2(1(1x0)0-1(0x1)0) 
((1x1)0 1 
((1x0)0 072 072 
((0x1)0 072 -0/2 

1(1 x 1) 2 ) 
((lxl) 2 | 1 



(C3) 
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The corresponding matrix elements ((S c q x S qS )s\(Scq x S 5 q)s) for open flavor decay are: 



|(0x0)o) |(lxl)o) 
((0x0)o| -1/2 V3/2 
((1 x 1) | -V3/2 -1/2 



(C4) 



l(lxl)i) 



V^72(|(l x 0)0 + |(0 x 1)0) ^l/2(\(l x 0)0 - |(0 x 1)0) 



((lxl)i| 
((1x0)0 
((0x1)0 



v/T72 








-1 



(C5) 



1(1 x 1) 2 ) 



(C6) 



((lxl) 2 | 



1 



In deriving the selection rules of Table IV, we bear in mind that the physical states for each J PC can contain 
any admixture of the various possible spin wave functions in Table III, and we thus consider only those selection 
rules which apply to all possible admixtures. 
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